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Abstract 

We present a linear time certifying algorithm that tests graphs for 3-edge- 
connectivity. If the input graph is not 3-edge-connected, the algorithm returns 
a 2-edge-cut. If the input graph is 3-edge-connected, the algorithm returns a con- 
struction sequence that constructs the graph starting from the graph with two nodes 
and three parallel edges using only operations that preserve 3-edge-connectivity. 
No previous algorithm returned a certificate in the case of a 3-edge-connected in- 
put graph. 

1 Introduction 

Many graph algorithms answer complex yes-no questions such as "Is this graph pla- 
nar?" or "Is this graph A:-vertex-connected?". These algorithms are not only nontrivial 
to implement, it is also difficult to test their implementations extensively, as usually 
only small test sets are available. It is hence possible that bugs persist unrecognized 
for a long time. An example is the implementation of the linear time planarity test of 
Hopcroft and Tarjan 0HT74I in LEDA IIMNU99I . A bug in the implementation was 
discovered only after two years of extensive use. 

Certifying algorithms BMMNSl II approach this problem by computing an addi- 
tional certificate that proves the correctness of the answer. This may, e.g., be either 
a 2-coloring or an odd cycle for testing bipartiteness, or a planar embedding or a Ku- 
ratowski subgraph for testing planarity. Certifying algorithms are designed such that 
checking the correctness of the certificate is substantially simpler than solving the orig- 
inal problem. Ideally, checking the coiTectness is so simple that the implementation 
of the checking routine (the checker) allows for a formal verification. In that case, the 
solution of every instance is correct by a formal proof IIABMRl II . 

Our main result is a linear time certifying algorithm for 3-edge-connectivity. We 
make use of a construction due to Mader that characterizes 3-edge-connected graphs 
and show how to compute this construction in linear time. If the input graph is not 3- 
edge-connected, a 2-edge-cut is computed. Previous algorithms 0GI9 1 1 |NI92| ITW092I 
ITsi07l|Tsl09J for deciding 3-edge-connectivity are not certifying; they deliver a 2-edge- 
cut for graphs that are not 3-edge-connected but no certificate in the yes-case. 

Our algorithm uses the concept of a chain decomposition of a graph introduced 
in BSchlOI . A chain decomposition is a special ear decomposition llLov85L They are 
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used in IISchl2bl as a common and simple framework for certifying 1- and 2-vertex, 
as well as 2-edge-connectivity. Further, |SchT2a| uses them for certifying 3-vertex- 
connectivity (an implementation is available in |Neul 1]). Chain decompositions are 
an example of path-based algorithms (see, e.g., Gabow IIGabOOII ). which use only the 
simple structure of certain paths in a DFS-tree to compute connectivity information 
about the graph. 

We use chain decompositions to certify 3 -edge-connectivity in linear time. Thus, 
chain decompositions form a common framework for certifying fc-vertex- and A;-edge- 
connectivity for k < 3 in linear time. We use many techniques from ||Schl2al . but 
in a simpler form. Hence our paper may also be used as a gentle introduction to the 
3-vertex-connectivity algorithm. 

2 Related Work 

Deciding 3 -edge-connectivity is a well researched problem, with many different linear 
time solutions known [019 1 NI92 TW092 Tsi07 Tsi09 1 . None of them is certifying. 

The paper IMMNSl II is a recent survey on certifying algorithms. For a linear 
time certifying algorithm for 3-vertex-connectivity, see |Schl2a, Neul 11. For general 
k, there is a certifying algorithm for A:-vertex connectivity in |LLW88| with running 
time 0{n^'^ +nk^-^) and error probability at most 1/n. There is a non-certifying algo- 
rithm OKarOOl for deciding A;-edge-connectivity in time 0{m\o^ n) with high probabil- 
ity. 

In IIGI91I . a graph transformation is described that transforms a graph G into a 
graph G' such that G is 3 -edge-connected if and only if G' is 3-vertex-connected. How- 
ever, the certifying 3-vertex-connectivity algorithm from IISchl2al INeul II is much 
more complex than the algorithm given here. Moreover, we were unable to find an 
elegant method for transforming the certificate obtained for the 3-vertex-connectivity 
of G' into a certificate for 3 -edge-connectivity of G. 

3 Preliminaries 

We consider finite undirected simple graphs G with n = \V{G) \ vertices and m — \E{G) \ 
edges and use standard graph-theoretic terminology from 0BMO8I . unless stated other- 
wise. 

A set of edges that leaves a disconnected graph upon deletion is called edge cut. 
For k> 1, let a graph G be k-edge-connected if n> k and there is no edge cut X in G 
with \X\ < k. Vertex connectivity is defined analogously. 

Let V -^G w denote a path P from vertex v to vertex w in G and let s{P) — v and 
t{P) = w be the source and target vertex of P, respectively (as G is undirected, the 
direction of P is given by s{P) and t{P)). Every vertex in f \ {s{P),t{P)} is called an 
inner vertex of P. 

Let T be an undirected tree rooted at vertex r. For two vertices x and y in T, let x 
be an ancestor of y and y be a descendant of xif x^V{r y). If additionally x^y, 
xisa. proper ancestor and y is a proper descendant. We write x <yifxisB.n ancestor 
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Figure 1: Construction of a 3-edge-connected graph using Mader-paths. Real vertices 
are depicted in black. The black edges exist akeady, while dotted grey vertices and 
edges do not exist yet. 

of y. The parent p{v) of a vertex v is its immediate ancestor The parent function is 
undefined for r. Let /T™ be the graph on 2 vertices that contains exactly m parallel edges 
and no self-loops. A decomposition of G is a set of subgraphs of G whose edge-sets 
partition E{G). 

Let subdividing an edge uv of a graph G be the operation that replaces uv with a 
path uzv, where z was not previously in G. Then a subdivision G' of a graph G is a 
graph obtained by a series of subdivisions. We call paths in G' that correspond to edges 
in G links. The real vertices of a subdivision are the vertices with degree at least three. 

In IIMad78l . Mader shows that all 3-edge-connected graphs can be constructed us- 
ing a small set of operations starting on a ■ 

Theorem 1 (Mader ||Mad78l ). Any 3-edge-connected graph can be constructed from 
a K2 using the following three operations: 

• Adding an edge (possibly parallel or a loop). 

• Subdividing an edge xy and connecting the new vertex to any existing vertex. 

• Subdividing two edges wx, yz and connecting the two new vertices. 

When we talk about subdivisions of 3-edge-connected graphs, the analogous op- 
erations add links. Edges are no longer subdivided as in the last two operations of 
Theorem [T] instead, links are added that end at inner vertices of other links, i.e. the 
vertex that would have been created by the subdivision is an inner vertex of a path that 
was added previously. We call the paths in this process Mader-paths. For an example 
construction, see Figure [T] 

Proposition 1. Let G be a subdivision of a 3-edge-connected graph and let C be a 
simple path disjoint from G except for its endpoints. If the endpoints of C are not both 
interior points of the same link of G, then C is a Mader-path for G. 

Proof. The only forbidden Mader-operation is subdividing the same link twice. □ 

In this paper we show how to find a Mader construction sequence for a 3-edge- 
connected graph in linear time. 

4 Chain Decompositions 

We use a very simple decomposition of graphs into cycles and paths. The decompo- 
sition unifies existing linear-time tests on 2-vertex- and 2-edge-connectivity ||Schl2bl 
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Figure 2: A chain C (solid black) with the green vertex being s{C), the red vertex being 
t{C) and the dashed chain being p{C). 

and is also the base structure for certifying 3 -vertex-connectivity efficiently. In this 
paper we show that it can also be used to find a construction sequence for a 3-edge 
connected graph. We define the decomposition algorithmically; a similar procedure 
that serves for the computation of low-points can be found in |Ram93 1. 

Let G be a simple but not necessarily connected graph and let T be a depth-first 
search forest of G. The DFS assigns a depth-first index (DPI) to every vertex. For 
every backedge e, s{e) and t{e) are the two end vertices of e such that s{e) is a proper 
ancestor of t{e) in T . 

We decompose G into a set = {Ci, . . . jCi-g-i} of cycles and paths, called chains, 
by applying the following procedure for each vertex v in ascending DFI-order: Let 
T' be the tree in the DFS -forest T that contains v and let r be the root of T' . For 
every backedge vw with s{vw) — v, we traverse the path w -^ji r until a vertex x is 
encountered that was visited before. The traversed subgraph vw U [w -^ji x) forms a 
new chain Ct with s{Ci) — v and f (C,) = x. Since every backedge defines one chain, 
there are precisely m — n + 1 chains. 

If X is not r, we set the parent of Q as p(C,) = Cp. We say a vertex v belongs to a 
chain C if vp{v) G C; by construction a vertex belongs to a unique chain. We define an 
order on the chains by setting C,- < Cj if / < j. See Figure|2] 

We call '^/a a chain decomposition. Clearly, C can be computed in time 0{n + m). 
We will use the following lemma. 

Lemma 1 ( ISchl2bll ). Let be a chain decomposition of a simple graph G. Then 

• G is 2-edge-connected if and only ifG is connected and the chains in partition 
E{G). 

• G is 2-connected if and only if it has minimum degree at least two and C\ is the 
only cycle in '£ . 
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5 The Algorithm 



An implementation in Python is publicly available at https : //github . com/ adrianN/edge- connectivity 
The overall structure of the algorithm is given in Algorithm [1] 

5.1 Computing the Certificate 



Algorithm 1 Certifying algorithm for 3 -edge connectivity 
procedure Connectivity(G=(V,E)) 

chains = MODIFIED_CHAIN_DECOMPOSITION(G) 

for chain in chains do 

for c in chain. children[type 2] do 

ADD(c) 

end for 

segments = [] 

for c in chain. children[type 1] do 
segments. append(Segment(c)) 
end for 

for c in chain. type3 do 

segment = MAKE_SEGMENT(c) 

if segment n chain. children[type 1] = then 

add_type3(c) 
end if 

end for 

ORDER_AND_ADD(segments) 
end for 
end procedure 



We show how to find a Mader sequence in linear time. The algorithm begins by 
computing a chain decomposition as described above and then modifying it slightly. In 
this phase we also check that the graph is 2-edge-connected using Lemma[T]and that it 
has minimum degree at least three. We will also use the following necessary condition 
for 3-edge-connectivity. 

Proposition 2. Let T be a DFS tree starting at r for a graph G. If G is 3-edge- 
connected, then the subtree of every child of r must be connected to r by at least two 
backedges. 

Using the chain decomposition, we can identify a subdivision in the graph as fol- 
lows. The first chain, C\, of the chain decomposition forms a cycle. By Proposition]!] 
there must be a chain that starts at r and has C\ as a parent, otherwise the graph is not 
3-edge connected. W.l.o.g. we can assume that C2 is this chain and f (C2) = x. We mod- 
ify the first chain of the chain decomposition by setting Q = x — s-j- r and C\=r — s^Cj x. 
Then Q U Ci U C2 form an initial K2 subdivision. Note that we only do this transforma- 
tion once, later chains might be cycles if the graph is not 2-vertex-connected. 
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Note that, except for r, every vertex belongs to the first chain that contains it. Start- 
ing with G3 = Co U Ci U C2, our goal is now to find either an order K on the remaining 
chains such that 

Gk= U C^i^j) 

0<j<k 

and C;t(i^) is a Mader-path for for every k G {3 , . . . , m — « + 1 } or to find a separating 
pair of edges for G if no such order exists. As the G, are constructed using Mader-paths, 
Theorem[T]implies that they are subdivisions of 3-edge-connected graphs. We say that 
we add a chain C, in step k if it is a Mader-path with respect to G^ and 7z{k) = i, and 
that a chain Q is part of the current graph in step k if C, C G^- 

We classify the chains {Ci, . . .C,„_„+i } into three types. Let C be a chain with 
parent C' = p{C). By construction, s{C) is a (not necessarily proper) descendant of 
s(C'). We classify C according to the location of s{C). There are three cases. 

• If s{C) belongs to C' , C has Type 1. Observe that t{C) —t't s{C) is contained in 
C' for type 1 chains. 

• If s{C) — s{C'), the chain is of Type 2. We also set the type of the chains starting 
at r and having parent Co to two, even though s{Cq) is not r. 

• Otherwise, i.e. if s{C) is an ancestor of f(C') and a proper descendant of i(C'), 
C is of Type 3. 

Let Gc be a union of chains. Gc is upwards-closed if, for any vertex v G Gc, the 
edge vp{v) is also part of Gc- 

Lemma 2. Let Gc be a upwards-closed union of chains that contains Co, Ci, and C2. 
Then any type two or three child of a chain in Gc is a Mader-path with respect to Gc- 

Proof. Let C be a chain in Gc- If C = Co, s{C) is real as Ci and C2 are in Gc- Otherwise, 
by the upwards-closedness of Gc and the construction of chains, s{C) is real as the edge 
s{C)p{s{C)) is in Gc and belongs to a chain different from C. Further, all children of C 
are ears with respect to Gc, that is, they are disjoint from Gc save for their endpoints. 

We apply Proposition[T] Let C' be a child of C. If C' has type two, s{C') = s{C) and 
hence C' is a Mader-path. If C' has type three, we have s{C) < s{C') < t{C) < t{C'). 
Since t{C) is real, t{C') — Sr s{C') is not a link. As Gc is a union of chains, no other 
path ? (C') -^Gc ''{C) can be a link either Hence C' is a Mader-path. □ 

We operate in phases. In phase /, / G [0, we first add all type two children of C/ 
and then all chains contained in segments with respect to C,. The segments with respect 
to Ci are sets of descendants of C, . They are defined by the following rules. 

• Every type one child of C; belongs to some segment. 

• Every type three chain C with s{C) G C, belongs to some segment. 

• If C belongs to some segment and p{C) does not belong to the current graph (i.e. 
the graph after adding all type two children of C,) then p{C) belongs to the same 
segment with respect to C, . 
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• That is all, i.e., only chains forced by one of the rules belong to a segment with 
respect to C, and segments are disjoint except if merged by the third rule. 

If these chains cannot be added, we will exhibit a separation pair 

As we shall see in Lemma [3] each segment is a subtree of the chain tree that con- 
tains no chains of type one, except possibly a type one child of C, as its root. In any 
segment the chains belonging to it are added parent-first, i.e., a parent is added before 
its children. In the sequel, we give details, in particular, we discuss in what order the 
segments are added. We refer to C, as the current chain of phase /. 

Proposition 3. The current chain is part of the current graph at the beginning of any 
phase. 

Proof. The initial current graph consists of chains Co,Ci and C2. If C; is of type one 
or two, it was added in the phase for its parent. If it is of type three, let s{Ci) lie on Cj. 
Then C, is in a segment of Cy, since j < i, Cj was added in the phase for Cj. □ 

Proposition 4. The current graph is always upwards-closed. 

Proof. Directly from the fact that for all chains that we add, the parent is already added, 
and the fact that we add nothing but whole chains. □ 

We now come to the details of a phase. By Lemma|2] type two children of C, can be 
added. We do so. Let Gc be the current graph after the addition of the type two chains. 
We will next discuss how to add the segments with respect to C, . Segments correspond 
to subtrees of the chain tree. Only the roots of the subtrees can be of type one as we 
next show. 

For a type three chain Do that is not part of the current graph and satisfies s{Do ) € Q, 
let the subsegment corresponding to Dq be all ancestors D/. < ■ • • < Dq of Dq that are 
contained in a segment of C,. 

Lemma 3 ( llSchl 11 Lemma 81]). LetD be a chain of type three with s{D) € Ci and let 
D' be any ancestor ofD such that p{D') is not part of the current graph. Then D' has 
type two or three. 

Proof. Assume otherwise, i.e. D' has type one. Then D' ^ D, since D has type three. 
Moreover, s{D') is a proper ancestor of s{D) as otherwise D would be of type two. 
Since D" = p{D') is not part of the current graph Gc, t{D') cannot be part of the current 
graph and hence, by upwards-closedness, t{D') must be a proper descendant of s{D). 
It follows that s{D) is an inner vertex of t{D') — ^j- s{D'). 

Since D' is of type one, D" contains the path t{D') -^j ^{D') and hence s{D) e D". 
Thus s{D) ^ Gc, a contradiction. □ 

In particular, if D' is of type one, it is the minimal chain in the segment and a child 
ofC,-. 

Corollary 1. All chains in a segment S can be added in parent-first order if its minimal 
chain can be added. 
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1 2 3 4 1 2 3 4 

Figure 3: Intervals for the blue segment with attachment points 1,2,4. 

Proof. By Lemma [3] all but the minimal chain in a segment are of type two or three. 
Hence the claim follows from Lemma|2] □ 

Next we describe how we compute segments. We only store segments that contain 
a type one child of C,, the others can be added immediately due to Lemma [2] and 
Corollary [U We compute the segments by traversing the subsegments of all type three 
chains C with s{C) G C,, starting with the chain having a source vertex with minimal 
DFI (i.e. the minimal chain w.r.t. <). If the minimal chain in the subsegment is not of 
type one, we immediately add all chains in the subsegment. Otherwise, we mark all its 
chains with the minimal chain of the segment. To achieve linear runtime, we stop the 
construction of a subsegment once a marked chain is reached. 

It remains to compute a proper ordering of the segments containing a type one 
chain or to exhibit a separation pair. For simplicity, we will say 'segment' instead of 
'segment containing a type one chain' from now on. 

For a segment S let the attachment points of S be all vertices in S that lie on C,. 
Note that the attachment points must necessarily be endpoints of chains in S and hence 
adding the chains of S makes the attachment points real. Type one children C of Q can 
be added if there are real vertices on t{C) -^j ■s(C'), therefore adding a segment can 
make it possible to add further segments. We reduce the problem of finding an order 
on the segments to a problem on intervals. 

W.l.o.g. assume that the vertices of C,- are numbered consecutively from 1 to \Ci\. 
For every segment 5, let aq < a i < . . . < a^. be the set of attachment points of S. We 
associate the following intervals with S 

{[aQ,a^]\l<t<k}yj{[ae,a,,]\l<l<k}, 

and for every real vertex v on C, we define an interval [0, v] . See Figure[3]for an example. 

We say two intervals [a, a'], [b,b'] overlap if a <b <a' <b' . Note that overlapping 
is different from intersecting; an interval does not overlap intervals in which it is prop- 
erly contained or which it properly contains. The overlap relation naturally defines an 
overlap graph for a set of intervals. For a set of segments Si,. . . ,Sk, let the reduced 
overlap graph be the graph on the segments and a special vertex for the real vertices 
in which there is an edge between two vertices if any intervals associated with these 
vertices overlap. 

Lemma 4. If the reduced overlap graph H for Ci is connected, we can add all segments 
ofCi. 
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Proof. Let R,Si,...,Sj be the vertices of i/ in a preorder, e.g. the order they are ex- 
plored by a DFS, starting at R, the vertex corresponding to the real vertices on C,-. An 
easy inductive argument shows that we can add all segments in this order Let C be 
the minimal chain of ^i. An edge between R and implies a real vertex on the path 
f(5'i) -^{Si), hence can be added by Corollary[T] Since the attachment points of 
a segment become real when it is added, a later segment 5, in the preorder can be added 
due to similar reasoning after its predecessors are added. 

If H is connected, this adds all segments. On the other hand suppose H is not 
connected and let ^ be a connected component that does not contain R. Let a be the 
minimal attachment point of any chain in JP^ and let b be the maximal one. Note that a 
and b are both different from i(C,) and f(C,), as these are real. 

Observe that all edges incident to vertices on b a belong either to C, or to a 
chain in . They can't belong to ancestors of Q, otherwise they would create real 
vertices, and they can't belong to a segment outside the connected component as they 
would induce an overlapping interval otherwise. 

As a and b are the extremal attachment points of y, the tree edge from a to its 
parent and the edge from Zj's predecessor on s{Ci) -^q t{Ci) to b disconnect a and b 
from the rest of the graph. □ 

It remains to show that we can find an order as required in LemmalU or a separating 
edge pair, in linear time. 

A naive approach that constructs the reduced overlap graph and runs a DFS will 
fail, since the overlap graph can have a quadratic number of edges. However, using a 
method developed by Olariu and Zomaya IIOZ96L we can compute a spanning forest 
of the graph in time linear in the number of intervals. The presentation in fOZ96 1 is for 
the PRAM and thus needlessly complicated for our purposes. A simpler explanation 
can be found in the appendix. 

Since the number of intervals created for a chain C, is bounded by 

|Childreni(Q)| +2|Type3(C0| + IVreailQ)!, 

where Childreni (C,) are the type one children of C,, Type3(C,) are the type three chains 
that start on Q, and VreailC, ) is the set of real vertices on C;, the total time spent on this 
procedure for all chains is 0{m). 

From the above discussion follows: 

Theorem 2. A Mader construction sequence can be found in time 0(n + m). 

5.2 Verifying the Certificate 

The certificate is either a separating edge pair, or a sequence of Mader-paths. For a sep- 
arating pair, we simply remove the two edges and verify that G is no longer connected. 

Checking the Mader sequence is slightly more involved. We assume that each edge 
in a Mader-path is doubly linked to the corresponding edge in G. Then it is easy to 
check that the Mader-paths partition the edges of G. 

Let G' be a copy of G. We begin by checking that the minimum degree of G' is 
at least three. Now we remove the Mader-paths again, in reverse order, suppressing 
vertices of degree two as they occur This can create multiple edges and loops. 
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Let G'j be the multi-graph before we remove the i-th path f,. At this point, f, must 
have been reduced to a single edge as inner vertices of P,- must have been reduced to 
degree two if P,- is an ear for G^. After removing the edge corresponding to P,, it must 
not be the case that both endpoints are still adjacent but have degree two, as then P, 
subdivided the same link twice. 

Lastly we check that G'2 is a K2 subdivision. 

6 Conclusion 

We presented a certifying linear time algorithm for 3-edge-connectivity based on chain 
decompositions of graphs. It is simple enough to use in a classroom setting and 
can serve as a gentle introduction to the certifying 3 -vertex-connectivity algorithm 

of IISchl2all . We also provide an implementation in Python, available at https : //github. com/adriEinN/edge-conne< 

There remain some open problems. Foremost, our algorithm only computes one 
separation pair. Is it possible to compute the 3-edge-connected components easily? 

Mader's construction sequence is general enough to construct ^-edge-connected 
graphs for any k>3, and can thus be used in certifying algorithms for larger k. So far, 
though, it is unclear how to compute these more complicated construction sequences. 
We hope that the chain decomposition framework can be adapted to work in these cases 
too. 
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A Computing a Spanning Subgraph of an Overlap Graph 



We first assume that all endpoints are pairwise distinct. We will later show how to 
remove this assumption by perturbation. 

For every interval / = [a,b] define its set of left and right neighbors: 

L{I) = {/' = [a',b'];a' <a<b' <b}, 
R{I) = {/ = [a',b'];a <a' <b< b'}. 

If the set of left neighbors is nonempty, let the interval /' € L{I) with the rightmost right 
endpoint be the immediate left neighbor of /. Similarly, if the set of right neighbors is 
nonempty, the immediate right neighbor of / is the interval in R{I) with the leftmost 
left endpoint. 

Lemma 5. The graph G' formed by connecting each interval to its left and right neigh- 
bor ( if any) forms a spanning subgraph of the overlap graph G and has exactly the 

same connected components. 

Proof. Clearly, every edge of G' is also an edge of G. For the other direction, assume 
/ and /' are overlapping intervals that are not cormected in G' and for which the left 
endpoint of / is as small as possible. Then a < a' < b < b', where / = [a,b] and 
/' = [a',b']. Since /' S /?(/), but / and /' are not connected, / has an immediate right 
neighbor J ^ I' . The left endpoint of J must thus be smaller than a' and the right 
endpoint of J must be larger than b (since / overlaps with J). If the right endpoint is 
smaller than b' , J and /' overlap. By repeating this argument for J and /', we must 
reach an interval U = [c,d] containing /'. Thus 

a<c <a' <b <b' <d. 

Starting from /' and going to left neighbors, we obtain in the same fashion an interval 
U' = [c',d'] with 

c' < a <a' <b < d' < b' . 

We conclude that U' and U overlap, but are not connected in G'. Since the left endpoint 
of U' is to the left of the left endpoint of /, this contradicts the choice of / and /'. □ 

It is easy to determine all immediate right neighbors by a linear time sweep over 
all intervals. We sort the intervals in decreasing order of left endpoint and then sweep 
over the intervals starting with the interval with rightmost left endpoint. We maintain a 
stack S of intervals, initially empty. If /i = [ai , fci] , . . . ,4 = [ak, bi^] are the intervals on 
the stack with I\ being on the top of the stack, then oi < 02 < . . . < fl/i and b\ <b2< 
... <bk,I\ is the last interval processed, and is the immediate right neighbor of 
if It has right neighbors. If does not have right neighbors, a^+i > b^. Let / = [a, fo] be 
the next interval to be processed. Its immediate right neighbor is the topmost interval 
on the stack with bi> b (if any). Hence we pop intervals from the stack while b> bi 
and then connect / to the topmost interval if b > an, and push /. The determination of 
immediate left neighbors is symmetric. 

It remains to deal with intervals with equal endpoints. We do so by perturbation. It 
is easy to see that the following rules preserve the overlaps-relation and eUminate equal 
endpoints. 
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(1) if a left and a right endpoint are at the same coordinate, then the left endpoint 
precedes the right endpoint. 

(2) if two left endpoints are equal, the one belonging to the shorter interval is smaller. 

(3) if two right endpoints are equal, the one belonging to the shorter interval is larger 

(4) if two intervals are equal, one is slightly shifted to the right. 

In other words, the endpoints of an interval = [a,b] are replaced by {{a, — l,b — a,i) 
and {b, l,b — a,i)) and comparisons are lexicographic. The perturbation need not be 
made explicitly, it can be incorporated into the sorting order and the conditions under 
which edges are added, as described in Algorithm|2] 



Algorithm 2 Finding a spanning forest of a overlap graph 

procedure SP(/ = {[ao,flo], • • • , [a? 
stack = [ ] 

sort I lexicographically in descending order 
for [l,r] in 1 do 

while stack not empty and r > top(stack) right endpoint do 
pop(stack) 

end while 

if stack not empty and r > top(stack) left endpoint then 

connect [l,r], top(stack) 
end if 

push(stack, [l,r]) 
end for 
stack = [ ] 

sort I lexicographically in ascending order where the key for [/, r] is [r, I] 
for [l,r] in 1 do 

while stack not empty and Z<top(stack) left endpoint do 
pop(stack) 

end while 

if stack not empty and I < top(stack) right endpoint then 

connect [l,r], top(stack) 
end if 

push(stack, [l,r]) 
end for 
end procedure 
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